Computation of Grobner Bases for Two-Loop Propagator Type Integrals 
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The Grobner basis technique for calculating Feynman diagrams proposed in [1] is applied to the two-loop prop- 
agator type integrals with arbitrary masses and momentum. We describe the derivation of Grobner bases for all 
integrals with IPI topologies and present elements of the Grobner bases. 



1. INTRODUCTION 

The most efficient algorithms to evaluate Feyn- 
man diagrams exploit recursive methods based on 
integration by parts technique [2] or technique of 
generalized recurrence relations [3] , [4] . It is easy 
to derive hundreds of recurrence relations but it 
is by far not so easy to formulate optimal algo- 
rithm how to use them for the reduction of a given 
type of integrals to a minimal set of basis inte- 
grals. In [1] it was proposed to use the Grobner 
basis method [5] as a mathematical background 
to solve this problem. This proposal can be real- 
ized in different ways. For example, the system 
of recurrence relations can be rewritten as a sys- 
tem of partial differential equations (PDE). The 
set of relations needed for reduction of Feynman 
integrals with different powers of propagators to 
the set of basis integrals will be Grobner basis of 
this overdetermined system of differential equa- 
tions. Information about the minimal basis of 
integrals also will be contained in the Grobner 
basis. Another possibility will be to rewrite recur- 
rence relations in terms of operators shifting pow- 
ers of propagators. Each equation should be con- 
sidered as an operator. Operators shifting powers 
of propagators satisfy an Ore algebra' condition 
and therefore, to find minimal set of recurrence 
relations from our system of equations ammounts 
to computation of it's Grobner basis in Ore alge- 
bra. 

Our preliminary investigation reveals that both 
approaches have some merits and shortcomings. 



The first approach certainly will be preferable in 
case when integrals depend on many mass scales. 
For integrals depending on two - three mass scales 
the second approach can be more efficient. 

In this article we present the computation 
of the Grobner basis for the two-loop propaga- 
tor type integrals with arbitrary momentum and 
masses. The derivation is based on a representa- 
tion of recurrence relations as partial differential 
equations. For integrals with three, four and five 
lines we used notation Jg , and F'^ respectively 
which was adopted in our paper [4]. 



2. TWO-LOOP 'SUNSET' INTEGRALS 

We begin our consideration with the so-called, 
'sunset' type of integrals given on Fig.l. 




Fig.l 

Recurrence relations for this type of integrals can 
be obtained from two equations [6], [2] 
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where 

ci = k'l - mf, C2 = (fci - fc2)^ - ml, 
C3 = {k2 - - ml, 
S^up = wigf_ii,qp + 

and Q^i^ is an arbitrary tensor. In our consid- 
eration we restrict ourselves to the second rank 
tensor in ki, k2 and external momentum q: 

Qfiv = xiki„ki^ + X2kijjk2p + X3ki^k2„ 

+Xik2vk2p + x^ki^q^ + xaki^qu (3) 
+X7k2i,qp + Xs,k2p,q„ + xgq^q^. 

Parameters Xi,Wj are arbitrary constants. After 
differentiation integrals with irreducible numera- 
tors 

J J C^C2Cz 

were expressed in terms of integrals with shifted 
dimension. For example, 
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Propagators with 'shifted' indices should be rep- 
resented as: 



1 d"" 1 
r! cJzJ Cj 



(5) 



where 



Oi0 give 28 different relations connecting 'sun- 
set' integrals with different shifts of the parameter 
of space-time dimension d and different products 
of one-loop tadpole type integrals. 

To find the Grobner basis of our system of PDE 
we used package Rif [7] which is distributed as 
part of Maple. 

Just for illustration we give a concrete example 
of the program written in Maple for computation 
of the Grobner basis. 

read 'equations': 
eq29:=zl*diff (Tl(zl) ,zl) 

- (d/2-l)*Tl(zl)=0: 
eqSO : =z2*dif f (T2 (z2) , z2) 

- (d/2-l)*T2(z2)=0: 
eqSl : =z3*dif f (T3 (z3) , z3) 

- (d/2-l)*T3(z3)=0: 
syst : ={seq(eq| I j , j=l. .31)}: 

syst : = syst union {eq29 , eq30 , eq31} : 
syst : = convert (syst , list) : 

with(DEtools) : 

invr := [zl,z2,z3]: 

dvar s : = [sun3 , sun2 , sunl , sunO , Tl , T2 , T3] : 
Rnk : = [ [0 , , , 150 , 100 , 50 , , , , 0] , 

[1,1,1, 0, 0, 0,0,0,0,0]] : 
bas := rifsimpC 

syst , dvar s , indep=invr , ranking=Rnk) ; 
GBasis:= bas [Solved]; 
quit ; 

First, the system of 28 differential equations 
stored in the file 'equations' is read in and then we 
add to it three differential relations for one-loop 
tadpole integrals with different masses: 



'Sunset' type integrals J3 with different shifts of 
the parameter of the space-time dimension d were 
considered as different functions, i.e: 



J? = sunO, 



= sunl. 



4+"^ = sun2, . 



(6) 



These substitutions allows one to transform the 
system of recurrence relations into a linear sys- 
tem of PDEs for the vector function P = 
{sunQ,sun\, . . At i>i — 1^2 — 1^3 — ^ Eqs. 
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In order to get rid of integrals with shifts of d 
we adopted here a ranking giving weight propor- 
tional to the shift of space time dimension. Also 
derivatives of functions have weight higher than 
the function itself. 

It took 5 min on a 1.6 GHz PC to obtain the 
Grobner basis consisting of 19 differential rela- 
tions for J3 with different shifts of d and three 
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relations for Ti{zi). The left hand sides of these 
19 relations are: 



92 J| 



52 J| 



92 J| 



52 J| 



dzj^' dzidz2' dzjdzg' dz2^ 

9V| 

dz2dz3 ' dzs"^ ' 



,2 7-<i+2 



52 J, 



3 



,2 Td+2 



52J|^ 
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li should be noted that keeping in Q only linear 
terms in /ci, fc2 gives incomplete basis for J3 . 

By using the Grobner basis reduction of 
J3 """^ . . . integrals with vi + V2 + > ^ can be 
done with elimination commands built in Maple. 
We wrote our own procedure for such a reduc- 
tion taking into account specific properties of our 
Grobner basis. This procedure works essentially 
faster than built in tools of Maple. 

3. INTEGRALS 



dzi'^dzs'^ 
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dzidz2dz3^ 



dz3'^dz2'^dz3 ' dzi '^dz2dz3 ' 



dzi^dz2dz3^' dzidz2^dz3^ ^'^^ 

In complete agreement with [4] the number of 
equations for J3 without shifts of d is enough to 
reduce integrals with vi + V2 + > 4 to four ba- 
sis integrals ,dJ^ / dzi,dJ^ / dz2,dJ^ / dz^, and 
three products of one-loop tadpoles. The basis 
for integrals with shifted d is not complete. As 
was shown in [4] these integrals can be expressed 
in terms of Jg and its first derivatives. Such a 
relation was not obtained from the 28 equations 
produced by Eqs. Q-ljSl)- ft turns out that if we 
shift (i^d-|-2in(Pl,(|21l and add the obtained 28 
equations to original 28 equations then from this 
enlarged system we will obtain in the Grobner ba- 
sis the required expression for J^'^^ which agrees 
with [4]. Therefore, computing Grobner basis we 
must include these equations in our system from 
the very beginning. 

Taking into account the relation for J^^^ in 
terms of J^,dJ^ / dzi,dJ^ / dz2,dJ^ / dz^, we will 
find that the Grobner basis for 'sunset' integrals 
consists out of 10 relations: 6 relations with sec- 
ond derivatives for Jg , one relation for jf^"^ and 
three relations for one loop tadpole integrals. Ex- 
pressions for derivatives of Jg"'"^, J^'^^, Jg"*"^ wiU 
be excluded from Q. 



In the same manner as it was done for the 'sunset' 
integrals we repeated computation of the Grobner 
basis for propagator type integrals with the topol- 
ogy given in Fig. 2. 




Fig. 2. 

The Grobner basis consists of ten relations which 
were already obtained for 'sunset' integrals, five 
relations for and three relations for 1-loop 
propagator integrals 
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New relations needed in the Grobner basis are 
dV^ dV^ dV^ dV^ 



(8) 
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As it was for the 'sunset' case the number of el- 
ements in the Grobner basis and the number of 
basic integrals for V type integrals is in agreement 
with [4]. 

4. INTEGRALS F"^ 

Grobner basis for F'^ integrals with the topology 
given in Fig. 3 includes basis of 'sunset' type in- 
tegrals, V type integrals and also new relations. 
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Fig.3. 

Six new relations for the following quantities 

Qpd Qpd Qpd Qpd Qpd ^^^^ 

dzi ' dz2 ' dz3 ' dz^ ' dz^ 

appeared in the Grobner basis of F'' integrals. 
Again these relations and the number of master 
integrals agree with the results given in [4] . 

5. CONCLUDING REMARKS 

The considered examples demonstrate that 
with already existing software one can compute 
Grobner bases for a rather complicated type of 
Feynman integrals. Our preliminary investiga- 
tion reveals also that Grobner bases for vertex 
integrals with arbitrary masses and external mo- 
menta can be computed with the existing com- 
puters and available software. Corresponding ex- 
pressions are rather long, sometimes thousands 
lines of codes. However after the Grobner basis 
is computed, it can be stored and used in differ- 
ent applications. At the present time we cannot 
use Grobner bases for reduction of integrals with 
arbitrary kinematics. The reason is that some 
masses or Gram determinants appear in denomi- 
nators and therefore one cannot make reductions 
when these factors are zero. Different solutions 
of this problem are possible and they will be pre- 
sented in future publications. 
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